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FRACTAL PROGRAMMINGES 


Cash from chaos 


Nick Brooks voyages into the strange and beautiful universe of the 
Mandelbrot Set — the mathematical expression of a revolutionary 
way of viewing nature 


Fractal geometry 


“Why is geometry often described as “cold” and “dry”? One reason 
lies in its inability to describe the shape of a cloud, a mountain, a 
coastline, or a tree. Clouds are not in spheres, mountains are not 
cones, coastlines are not circles and bark is not smooth, nor does 

lightning travel in a straight line. More generally I claim that many 
patterns of nature are so irregular and fragmented that, compared 
with Euclid — a term used in this work to denote all standard Ge- 
ometry — Nature exhibits not simply a higher degree, but an alto- 

gether different level of complexity. The number of distinct levels of 
natural patterns being for all practical purposes infinite”. 
Benoit B Mandelbrot — The Fractal Geometry of Nature 


At the heart of chaos theory lies the Mandelbrot Set — an 
immensely complex mathematical shape which expresses a whole 
new way of understanding natural phenomena 


he Mandelbrot Set lies 
like a black hole, a vortex 
of infinite complexity at 
the centre of an imagi- 
nary two-dimensional sea of 
numbers called the complex plane. 
From a distance, the Set could be 
a psychoanalyst’s inkblot, sur- 
rounded by multicoloured waves 
lapping against its shores. 

It seems strange that this bio- 
logical-looking shape should lie at 
the heart of a new science, and a 
new attitude to the way we under- 
stand and predict natural phenom- 
ena — chaos theory. But this is the 
one and only Mandelbrot Set — the 
most complex mathematical ob- 
ject ever seen. 

There is no mysterious alchemy 
at work here. Each point on this 
complex plane has been tested 
against a terse algorithm many 


tens, hundreds or thousands of 
times. Inside the Set, the points 
which have been captured and are 
unable to escape after repeated 
testing are displayed as black. 
Outside the Set are the points 
which have escaped before testing 
was completed, having near infi- 
nite values, and they are displayed 
as coloured. 

The Set is named after Benoit 
B Mandelbrot, a mathematician 
and research fellow from IBM’s 
Research Centre at Yorktown 
Heights, New York. Dating from 
1979, the Set has risen to become 
something of an emblem of chaos 
and the brightest star in modern 
mathematics. It was Mandelbrot’s 
explorations into the world of 
geometry which led to this and 
other discoveries, which he named, 
collectively, fractal geometry. 


How long is the coastline of 
Britain, and is it even possible to 
accurately measure it? Think about 
it — if you were to examine the 
southern coast of Cornwall from a 
satellite, and then examine a part 
of the same coastline on your 
knees, would one seem less com- 
plex than the other? No. The pat- 
tern of bays and peninsulas visible 
from space would be present in a 
microscopic scale as you leant over 
the side of a rock pool. 

If you were then to examine a 
square millimetre of the edge of the 
rock pool, you would again see a 
complexity of pattern and form 
equalled by the view from space. 
To all intents and purposes, this 
macro-to-micro scale of complex- 
ity of pattern and form is infinite. 
The closer you examine a natural 
object, the more detail is visible. It 
is the existence of this seemingly 
unmeasurable and recursive self- 
similarity in nature that we see at 
the boundary or shore of the 
Mandelbrot Set. 


Micro-worlds 

With the aid of a simple program, 
we convert our PC into a micro- 
scope and start our voyage into the 
complex micro-world of the Man- 
delbrot Set. Our voyage takes us 
closer to the boundary separating 
the Set from the sea of colours to 
reveal an intricate and minute 
coastline of instability. Shapes 
looking like seahorses’ tails and 
islands of molecules appear in 
microscopic detail. Shooting fur- 
ther into the crook of one of these 
tails, we see new patterns like solar 
flares and frost crystals loom into 


being from seemingly chaotic ar- 
eas of colour. We travel further, 
aiming towards an area of black 
embedded like a dark jewel in a 
spiral of colour. Approaching it, 
we see the area is a sibling set of 
the Mandelbrot Set, many thou- 
sands of times smaller, but still 
familiar. It is equally as beautiful 
and complex, but never identical, 
with its own coastline of instabil- 
ity and chaos. 

This is an infinite world. There 
is no end to the recursion of com- 
plexity of pattern, form and self- 
similarity on these alien 
coastlines.We will never find a 
repetition or a straight line, which 
— if it existed — would mark the 
end of our fractal journey. 


Programming 

The Mandelbrot program needs 
just a few essential pieces. The 
central engine of the program is a 
remarkably simple algorithm: 


z->z7+Cc 


where z is a complex number that 
is allowed to vary and c is a fixed 
complex number. If, like me, you 
hold about as much love for alge- 
bra as you hold for the Commu- 
nity Charge, don’t get too uncom- 
fortable — this is easy! The algo- 
rithm is a simple loop which takes 
a value for z, then squares it, adds 
c and then passes the result back 
to z, to give a new starting value 
to be looped again and again — 
iterated, in other words. 


Complex numbers 
At this point we should take a 
quick look at what is meant by 
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JOURNEY INTO INFINITY 


Just as the patterns visible in a coastline seen from space are reproduced in a microscopic portion of the edge of a rock pool, so the 
Mandelbrot Set reproduces itself infinitely. Figures A to F show progressive magnifications of portions of the Set until, in the middle of 
figure F, we can see a tiny sibling set of the Mandelbrot Set itself emerging 
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THE MANDELBROT SET 


complex numbers. Created and 
agonized over by 19th Century 
mathematicians, complex numbers 
were created to combine real 
numbers with an imaginary value 
of sqrt.—1. Interestingly, they area 
very good example of how pure 
and abstract 19th Century mathe- 
matical theory was, and have found 

“~~ remarkable applications in today’s 
electronic technology. 

A complex number is made of 
two parts which, for purely his- 
torical reasons, are called real and 
imaginary. For example, 2+3i is a 
complex number with real part 2 
and imaginary part 3i. The i next 
to the 3 has a value of sqrt.—1, and 
also shows which part of the 
number is imaginary. The modern 
view of complex numbers is that 
they are most simply defined as 
points or co-ordinates in a plane — 
the complex plane. Thus 2+3i is 
the address for a point 2 east and 
3 north on the complex plane. 

Adding complex numbers (ad- 
dition and multiplication are used 
in the algorithm) is simple. Just 
add the real parts together and the 
imaginary parts together: 


um 2+4i 
+ 9-2i 
= 11+2i 


Multiplication is as simple, re- 
membering that i squared = —1: 


2+3i 
xX 24+3i 
6i+9i2 
4+6i 
4+12i+9i7 
4+12i-9 
= -54+12i 


When complex numbers are ap- 
plied to z —> z* + c iteratively, 
strange things begin to happen to 
particular values of c. For ex- 
ample, here is what happens when 
Cis A+: 

First iteration: 1+3i 
Second iteration: -7+7i 
Third iteration: 1+97i 


Note that the real and imaginary 
parts may grow, shrink or change 
sign after each iteration. If this 


process continues, the resulting 
values for z will get progressively 
larger, and very shortly will scream 
off into infinity. 


Border skirmishes 

Obviously, when programming the 
Set, we must have some means of 
breaking out of the z —> 2 +c 
loop. Otherwise, after just a few 
iterations, we will exceed the ca- 
pacity of any computer with near 
infinite values for z. So we adopt 
two methods for breaking out of 
the loop, or bounding. 

1) The program must watch its 
running total. If the value for z 
exceeds 2, the theory of complex 
number iteration guarantees us 
that further iteration will drive z 
off into infinity — and therefore z 
can be assumed to be outside the 
Set and the program will move to 
the next point to be tested. 

2) We can simply limit the 
amount of iterations the loop can 
make, thereby saying that any 
point on the complex plane which 
is still <= to 2 after this iteration 
limit can be assumed to be a 
member of the Set, and again the 
program will move on. 

After each point on the com- 
plex plane has been tested accord- 
ing to these simple rules, the pro- 
gram displays the result graphi- 
cally. All points inside the Set are 
coloured black. Those points out- 
side the Set are assigned colours 
according to how soon they break 
out of the loop. For example, 
those points breaking out after 
three iterations might be displayed 
as green, after four red, after five 
blue, and so on. A zoom lens for 
the Mandelbrot is the next essen- 
tial thing any would-be fractal 
traveller should create for the 
program, without which the true 
beauty and intricacy of the Set will 
remain hidden. This a relatively 
simple task. Generate two cursors 
which read values for the ranges 
and then input these as the para- 
meters for a new drawing — creat- 
ing a fully interactive micro-view- 
port to the Set. The power of this 
zoom lens will depend to some 


(In Borland Turbo C) 


#include <graphics.h> 
#include <math.h> 
#include <conio.h> 
#include <alloc.h> 
#include <stdio.h> 
#include <dos.h> 


main() 


esscreen*/ 
for(i=O; i<=3; i++) 


cr=xrange+k*4/(getmaxx()/scale); 
ci=yrange+j*3/(getmaxy()/scale); 
zm=zr=Zi=zr2=zi2=0.0; 

/* Mandelbrot engine */ 
for(it=-1; (it<iters) && (zm<4); it++) 
{ 


z1=zr2-zi2+cr; 
22=2*zr*zi+ci; 


zr=21; 
zi=z2; 
zr2=zr*2r; 
zi2=zi*zi; 


if (it==iters) 


} 
} 
} 
} 


.S 
off 


degree on the effective size of 
numbers that your PC can ma- 
nipulate. An Amiga will produce 
magnifications in the order of 
x800,000. If the Amiga was a real 
microscope, its power would equal 
that of an electron microscope. 

Other useful additions to your 
program could include a screen- 
save procedure to save your frac- 
tal snapshots for easy retrieval 
without lengthy recalculation and 
redrawing. Or perhaps giving the 
colour values an additional height 
value which, when combined with 
a perspective viewing procedure, 
will create a contoured rock relief 
of the Mandelbrot Set. 

But the most compelling and 
aesthetically-pleasing images are 
created by careful manipulation of 


int graphdriver=DETECT, graphmode, iters=32, i, j, k, it; 
float scale = 1.0, xrange = -2.5, yrange = -1.5; 
float cr, ci, zr, Zi, zm, 21, 22, zr2, 2i2; 
int a[4]={0, 1, 0, 1}, b[4]={0, 1, 1, 0}; 


initgraph( &graphdriver, &graphmode, “”’); /*initialis- 


for(j=a[i]; j<=getmaxy()/scale; j+=2) 


for(k=b[i]; k<=getmaxx()/scale; k+=2) 
{ 


zm=zr2+zi2; 


/* drawing commands */ 


putpixel(k, j, (int)(fmod(it+2, 16))); 
if (kbhit()) exit(); 


the colour palette. Examples of 
these intricate images can be seen 
in the horde of articles and books 
published on fractals. »88 


RECOMMENDED READING 


The Science of Fractal 
Images — Peitgen and 
Saupe (Springer-Verlag). 
The Beauty of Fractals — 
Peitgen and Richter 


(Springer-Verlag). 
Chaos — James Gleick 
(Cardinal). 

The Fractal Geometry of 
Nature — B Mandelbrot 
(Freeman). 


Fractals are everywhere. They are 
all around us. In the texture of 
bark, in the shape of clouds and in 
the flickering of a flame. To a non- 
scientist, it may seem odd that 
such familiar things as these have 
become the focus for such intense 
research. But familiarity is not 
enough to ensure that scientists 
have the tool for an adequate 
understanding. In fact, many sci- 
entists and thinkers have become 
aware of the inadequacies of the 
established wisdom in areas as 
seemingly fundamental as natural 
geometry. Fractals answer many 
of these questions which, all too 
often, had been shied away from 
before because they were consid- 
ered too complex, and therefore 
unsolvable. 

The benefits of chaos theory 
are now being felt strongly in a 
complete sweep through many dis- 
ciplines and sciences. Fractals will 
soon be used to compress graphics 
data into small enough packages 


to pass along telephone lines, 
making possible new forms of 
visual telecommunication such as 
the videophone. In other areas of 
research, such as economics, fluid 
dynamics, meterology and anthro- 
pology, chaotic algorithms are 
being used to model natural be- 
haviour graphically, with an at- 
tention to detail not available 
before. They can tell us more 
about how well we can reliably 
predict, for example, economic 
trends, as opposed to building 
simplistic models which pay little 
or no attention to the chaotic 
behaviour present in many large 
and even small scale systems. 

In films, fractals are being used 
to generate special effects, as in 
Star Trek’s Genesis effect — the 
cosmic rebirth of a planet — which 
used fractal procedures to create 
the fire and mountains for the se- 
quence. 

It is safe to predict that most 
graphics software and hardware 
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“Art is a lie that lets us recognize the truth.” 
Pablo Picasso 


will soon have fractal procedural 
primitives for the manipulation of 
surfaces, enabling users to create 
natural-looking models of moun- 
tains, water, clouds, trees and so 
on, all in the quest for realism (the 
Symbolics Maclvory Graphics 
Workstation already has many of 
these facilities). 

The Mandelbrot Set is the most 
complex of all these fractals, and 
really is responsible for giving chaos 
a face. Just to hear a scientist 


admiring its frightening complex- 
ity and intricate beauty, more in 
the terms of Blake and Byron, 
forces us to question whether it is 
art or science — two things that we 
are traditionally told are as differ- - 
ent as intuition and analysis. The 
answer is simple. It is both. 


Benoit Mandelbrot always takes 
money to the pub with him 


SS 


4 COMPUTERS 


(All genuine UK models) 
COMMODORE 


Amiga B2000 + Colour Monitor ........... 
Amiga 2088 PC/XT Bridgeboard + 5.25” Drive 
Amiga 2286 PC/AT Bridgeboard + 5.25” Drive 


AMSTRAD 


Amstrad PC 1640 DD Mono .... 
Amstrad PC 1640 DD Colour 
Amstrad PC 1640 DD EGA... 
Amstrad PC 1640 HD30 Mono. 


YOAYGCAG 


Amstrad PC2086 DD Colour.......... 
Amstrad PC2086 DD 12” EGA/VGA 
Amstrad PC2086 HD30 12” EGA/VGA .. 
Please call for models not listed. 

ATARI 

Atari 520STFM Explorer Pack ............... 
Atari 520STFM Power Pack with 20 gai 
Atari 520STFM as above + Mono Monitor . 
Atari 1040STFM + Superbase/VIP Pro/MS 


SBC (inc On-Site Maint) 
SBC HD30 Single Drive + 30Mb Drive... 


SBC HDD40 Single Drive + 40Mb Drive . 
SBC FDC Twin Drive CGA Colour........... 


Please call for models not listed 


bABbABAEBAC ABC AG 


B2000 + 1084 + XT Bridgeboard + 20Mb Hard Drive 
B2000 + 1084 + At Bridgeboard + 20Mb Hard Drive 
B2000 + 1084 + At Bridgeboard + 40Mb Hard Drive 
Amiga B2000 only........6..........00005 


Amiga 500 Batman Pack Interceptor/Batman/DPaint II/NZ Story 
Amiga 500 Batman Pack + Philips 8833 Stereo... . 
Commodore 64C New Light Fantastic Pack... .. . . 


Amstrad PCW9512 inc Cut Sheet Feeder ....... . 


Amstrad PC2086 SD Mono 2086's Inc Lotus Symphony 


Atari 1040STFM Mono System + aboves/w...... 
SBC FD Twin Drive 640K 8MHz HDA/CGA Clock/Cal . 


SBC HD30C Single Drive + 32MB Drive CGA Colour . re 
SBC 286SMD MUltidrive 20Mb + 5.,25” + 3.5” 12 MHz Mono 1469.00 


HAND HELD 


AAMTANATANHAAANAANANAANAAANANATAANATATNA 
mY OFFICIALDEALER * AMIGA % AMSTRAD x ATARI 


MODEMS 
Psion Organiser CM/XP with 2-line display........... 65.17/86.91 V21/22/23/22bis + § 7 
1216.52 | Psion Organiser LZ/1.264 with 4-line display . 108.65/139.09 | Amstrad MC2400 internal PC V21/22/23/22bis + SW 7800 
trees . Psion MC200/400/600 (New Laphelds) ... vce PHONE Amstrad SM2400 External V21/22/23/22bis .............. . 229.00 
Poon 290-00: ‘psion HardWwaralSottware Add-ons: Pace Linnet Internal PC V21/23 AA AD + LED's + S/W 99,00 
Ce gored Psion Organiser 80 Col Printer Batt/Mains Pace Linnett External V21/23 AA AD . 130.39 
‘ . .PHONE 
. 216.52 | PRINTERS SAGE BUSINESS SOFTWARE 
anh ayaae Citizen 120D Parallel or Commodore .................... 113.00 | Mainlan Network Starter H/ware & S/ware for 2 PCs.. - 350.00 
51049 SFLOOPaaRA ftir ng b. pec unee eee. 138.26 Financial Controller Accounts: Sales/Pur/Inv/Stock .. 289.00 
pe ae 121.73 | StarLC-10 Colour Parallel . 178.26 once Plus .. “ bey 
pune ie i Star LC-10 Il Faster LC-10 ........secseeeee 160.87 ccountant.. . . 
Star LC24-10 24 pin Multi-font 170/57 cps .. . .. 207.83 Bookkeeper.. + 75.00 
beeeeeees 451.30 | Star FR10/15 inc on site maint (New) 250/63 cps..... . 319/419 | Payrollll...... - 109.00 
- 534.00 | Star XB24-10/15 inc. on-site main (New) 200/67 cps .. . 419/549 | PC Write Word Processor 75.00 
. 622.00 PC Planner...... 59.00 
f aes MONITORS PC Planner Il. ped 
° ‘ Retrieve Ill Relat - 109. 
. ae Atari High Res Mono Monitor... 2.0.0.6. ..0eecee ees 119.00 
: . Philips 8833 14” col RGB/Comp Med Res Stereo 207.83 
- 799.00 | Philips 7502 Green 12” Comp Video + Sound............. 85.22 | MISCELLANEOUS SOFTWARE 
979.00 Learn to: Type/Use Your PC/Use DOS, each 21.70 
DISK DRIVES SimpAcct low cost full function PC accounts... 89.00 
233.91 Amiga 2000 Internal 3.5” Disk Drive.................005 60.83 
312.17 20Mb Disk Drive for PC/A2000 inc Controller . . 156.52 MISCELLANEOUS 
- 390.00 20 Mb SCSI Disk Drive for A2000........ ; Amiga Pace Modem Package V21/23 Prestel, BBS etc inc s/ware 207.83 
- 390.00 Amiga 500 AMDRIVE 30Mb Disk Drive. Amiga 2000 8Mb RAM Board Populated 2Mb (Compatible)........ 346.96 
cece y 468.26 Amiga A590 20Mb Autoboot for A500 ae Amiga 512K RAM/Clock eres Ane i é 86.91 
FRAWA EXDONUMNG 525,532 secravsie tis ssasonsebiarsekan hates ketyacdsooncds 343.4 Disk Head Cleaner Kits — 3”/3.5"/5.25” .. - 6.08/5.64/5.21 
649.00 Microbotics Autobooting 32/48Mb Amiga Disk Drive .. maa ets MinIGEN Genlock for Amiga Bs ag 
. 779.00 Supra cOMb tor Atm Sh. heaves oe ort suan Tbek's +.. 399.1 Printer 2-way T-switch....... aa PODS 
. 869.00 Cumana 3.5” Disk Drive for A500 with Disable & Pass Thru........ 73.91 SBC 500 Watt Uninterruptible PSU « 359.00 
. 799.00 Cumana CSA354 1Mb for Atari ST : 82.61 SBC 1000 Watt Uninterruptible PSU . + 799.00 
+ 1039.00 Commodore 1541-11 for C64/128... v= 120,87 Surge Protector 13APlug............ aie hleee 
Scribecard 40Mb PC HardCard.............. 0.0 c ee seee 249.00 Surge Protector 3-Way Adaptor/4-Way Dist Unit ....... 17.35/13.87 


Please add 15% VAT to above prices. CARRIAGE £5 (Express £10). Software & Small Items free. Prices subject to change without notice. 


8 Ruswarp Lane, WHITBY, N. Yorks YO21 1ND 


Tel: 0947 600065 (9am-7 
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